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\ Abstract. Motivated by a quaternionic formulation of quantum mechanics, we discuss 

^— ^ , quaternionic and complex linear differential equations. We touch only a few aspects of 

^ ■ the mathematical theory, namely the resolution of the second order differential equations 

f*^ | with constant coefficients. We overcome the problems coming out from the loss of the 

fundamental theorem of the algebra for quaternions and propose a practical method to 

i-^ ■ solve quaternionic and complex linear second order differential equations with constant 

coefficients. The resolution of the complex linear Schrodinger equation, in presence of 

f-H , quaternionic potentials, represents an interesting application of the mathematical mate- 

* ■ rial discussed in this paper. 

03 

B 
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: 1 INTRODUCTION 

There is a substantial literature analyzing the possibility to discuss quantum systems by adopting 
quaternionic wave functions Jl|-p[. This research field has been attacked by a number of people 
leadin g to substantial progress. In the last years, many articles |l5|-|3l|, review papers ]32]-|34|| and 
books [p5|-^7| provided a detailed investigation of group theory, eigenvalue problem, scattering theory, 
relativistic wave equations, Lagrangian formalism and variational calculus within a quaternionic for- 
mulation of quantum mechanics and field theory. In this context, by observing that the formulation of 
physical problems in mathematical terms often requires the study of partial differential equations, we 
develop the necessary theory to solve quaternionic and complex linear differential equations. The main 
difficulty in carrying out the solution of quaternionic differential equations is obviously represented 
by the non commutative nature of the quaternionic field. The standard methods of resolution break 
down and, consequently, we need to modify the classical approach. It is not our purpose to develop 
a complete quaternionic theory of differential equations. This exceeds the scope of this paper. The 
main objective is to include what seemed to be most important for an introduction to this subject. 
In particular, we restrict ourselves to second order differential equations and give a practical method 
to solve such equations when quaternionic constant coefficients appear. 

Some of the results given in this paper can be obtained by translation into a complex formal- 
ism [To], [Il| . Nevertheless, many subtleties of quaternionic calculus are often lost by using the 
translation trick. See for example, the difference between quaternionic and complex geometry in quan- 
tum mechanics |p2|,p4[|, generalization of variational calculus |^,|l^], the choice of a one-dimensional 
quaternionic Lorentz group for special relativity [glj , the new definitions of transpose and determi- 
nant for quaternionic matrices j2Sj ] . A wholly quaternionic derivation of the general solution of second 
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order differential equations requires a detailed discussion of the fundamental theorem of algebra for 
quaternions, a revision of the resolution methods and a quaternionic generalization of the complex 
results. 

The study of quaternionic linear second order differential equations with constant coefficients 
is based on the explicit resolution of the characteristic quadratic equation |38-41|. We shall show 
that the loss of fundamental theorem of the algebra for quaternions does not represent a problem in 
solving quaternionic linear second order differential equations with constant coefficients. From there, 
we introduce more advanced concepts, like diagonalization and Jordan form for quaternionic and 
complex linear matrix operators, which are developed in detail in the recent literature |2^-|3l|| and we 
apply them to solve quaternionic and complex linear second order differential equations with constant 
coefficients. 

As application of the mathematical material presented in this paper, we discuss the complex linear 
Schrodinger equation in presence of quaternionic potentials and solve such an equation for stationary 
states and constant potentials. We also calculate the relation between the reflection and transmition 
coefficients for the step and square potential and give the quaternionic solution for bound states. 

This work was intended as an attempt at motivating the study of quaternionic and complex linear 
differential equations in view of their future applications within a quaternionic formulation of quantum 
mechanics. In particular, our future objective is to understand the role that such equations could play 
in developing non relativistic quaternionic quantum dynamics Q| and the meaning that quaternionic 
potentials psjp^fl could play in discussing CP violation in the kaon system J36|. 

In order to give a clear exposition and to facilitate access to the individual topics, the sections 
are rendered as self-contained as possible. In section |2| we review some of the standard concepts used 
in quaternionic quantum mechanics, i.e. state vector, probability interpretation, scalar product and 
left/right quaternionic operators [^9[^,^2 45 1. Section || contains a brief discussion of the momentum 
operator. In section ||, we summarize without proofs the relevant material on quaternionic eigenvalue 
equations from pH . Section ^ is devoted to the study of the one-dimensional Schrodinger equation 
in quaternionic quantum mechanics. Sections ^ and R provide a detailed exposition of quaternionic 
and complex linear differential equations. In section we apply the results of previous sections to 
the one-dimensional Schrodinger equation with quaternionic constant potentials. Our conclusions are 
drawn in the final section. 



2 STATES AND OPERATORS IN QUATERNIONIC QUANTUM MECHANICS 

In this section, we give a brief survey of the basic mathematical tools used in quaternionic quantum 
mechanics [^2^3^|. The quantum state of a particle is defined, at a given instant, by a quaternionic 
wave function interpreted as a probability amplitude given by 

*(r) = [/„ + h-/](r), (1) 

where h = (i, j, k), f = (f ± , / 2 , / 3 ) and / m : R 3 — > M, m = 0, 1, 2, 3. The probabilistic interpretation of 
this wave function requires that it belong to the Hilbert vector space of square-integrable functions. 
We shall denote by J- the set of wave functions composed of sufficiently regular functions of this vector 
space. The same function ^(r) can be represented by several distinct sets of components, each one 
corresponding to the choice of a particular basis. With each pair of elements of J- ', ^f(r) and $(r), we 
associate the quaternionic scalar product 

0,$) = J d 3 r *(r)$(r) , (2) 

where 

*(r) = [/„-fc-/](r) (3) 

represents the quaternionic conjugate of \t(r). 

A quaternionic linear operator, O b , associates with every wave function ^/(r) € T another wave 
function O a ^(r) S J-, the correspondence being linear from the right on H 



O a [■$ 1 {r)q 1 +y 2 (r)q 2 } = [O a *i(r)] «i + [0„tf a (r)] ^ , 
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<Zi,2 G H. Due to the non-commutative nature of the quaternionic field we need to introduce complex 
and real linear quaternionic operators, respectively denoted by O c and O m , the correspondence being 
linear from the right on C and R 

O c [*!(r)zi + * a (r)^] = [O c *i(r)] * + [O c tf a (r)] z 2 , 
O h [*,(r)\+* 2 (r)\,] = [0 H *i(r)] A t + [0 H * a (r)] A 2 , 

z L2 € C and A 1-2 G M. 

As a concrete illustration of these operators let us consider the case of a finite, say n-dimensional, 
quaternionic Hilbert space. The wave function W(r) will then be a column vector 



\*»7 



Quaternionic, complex and real linear operators will be represented bynxn quaternionic matrices 
M n [A O] , where O represents the space of real operators acting on the components of W and 
A = (A E , A c , A K ) denote the real algebras 

An : { 1 , L , i? , £*i?} 16 , 
A : {1 , i , Ri, LR>} 8 , 
A s : { 1 , £ } 4 , 



generated by the left and right operators 

L := (L i: Lj, L k ) , R := (R i: R k ) 

and by the mixed operators 

L * R := {L p R q } p,q = i,j,k. 
The action of these operators on the quaternionic wave function \& is given by 

= m , R^ = ^h . 
The operators L and R satisfy the left/right quaternionic algebra 

L\ — L 2 . = L k = L { LjL k = R 2 = R^ = R 2 k — R k RjRt = — 1 
and the following commutation relations 

[L p ,R q ]=0. 



(4) 
(5) 



3 SPACE TRANSLATIONS AND QUATERNIONIC MOMENTUM OPERATOR 

Space translation operators in quaternionic quantum mechanics are defined in the coordinate repre- 
sentation by the real linear anti-hermitian operator [ [36| 

d = (d^,d y ,d z ) . (6) 

To construct an observable momentum operator we must look for an hcrmitian operator that has 
all the properties of the momentum expected by analogy with the momentum operator in complex 
quantum mechanics. The choice of the quaternionic linear operator 

V L = -L % hd, (7) 

as hermitian momentum operator, would appear completely satisfactory, until we consider the trans- 
lation invariance for quaternionic Hamiltonians, 7i q . In fact, due to the presence of the left acting 
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imaginary unit i, the momentum operator (|7]) does not commute with the j'/fc-part of Ti q . Thus, 
although this definition of the momentum operator gives an hermitian operator, we must return to 
the anti-hermitian operator d to get a translation generator, [d, H q ] — 0. A second possibility to be 
considered is represented by the complex linear momentum operator, introduced by Rotelli in |^| , 

V R = -R i hd. (8) 

The commutator of V R with a quaternionic linear operator O b gives 

[v R ,o]* = h[0,d] m . 

Taking O u to be a translation invariant quaternionic Hamiltonian TL q , we have 

[V*,H t ] =0 . 

However, this second definition of the momentum operator has the following problem: the complex 
linear momentum operator T R does not represent a quaternionic hermitian operator. In fact, by 
computing the difference 

which should vanish for an hermitian operator "P R , we find 

(tf , Ph*) - (Vr*, *) = h [i, {% 0$)] , (9) 

which is in general non-vanishing. There is one important case in which the right-hand side of equa- 
tion (|]) does vanish. The operator "P R gives a satisfactory definition of the hermitian momentum 
operator when restricted to a complex geometry [ff5f , that is a complex projection of the quaternionic 
scalar product, ('J, "P R <I ) ) C . Note that the assumption of a complex projection of the quaternionic scalar 
product does not imply complex wave functions. The state of quaternionic quantum mechanics with 
complex geometry will be again described by vectors of a quaternionic Hilbert space. In quaternionic 
quantum mechanics with complex geometry observables can be represented by the quaternionic hermi- 
tian operator, H , obtained taking the spectral decomposition pl"|] of the corresponding anti-hermitian 
operator, A, or simply by the complex linear operator, — AR it obtained by multiplying A by the 
operator representing the right action of the imaginary unit i. These two possibilities represent equiv- 
alent choices in describing quaternionic observables within a quaternionic formulation of quantum 
mechanics based on complex geometry. In this scenario, the complex linear operator T R has all the 
expected properties of the momentum operator. It satisfies the standard commutation relations with 
the coordinates. It is a translation generator. Finally, it represents a quaternionic observable. A review 
of quaternionic and complexified quaternionic quantum mechanics by adopting a complex geometry 
is found in Q]. 



4 OBSERVABLES IN QUATERNIONIC QUANTUM MECHANICS 

In a recent paper fl3l| , we find a detailed discussion of eigenvalue equations within a quaternionic 
formulation of quantum mechanics with quaternionic and complex geometry. Quaternionic eigenvalue 
equations for quaternionic and complex linear operators require eigenvalues from the right. In partic- 
ular, without loss of generality, we can reduce the eigenvalue problem for quaternionic and complex 
linear anti-hermitian operators A G M n [A B ® O] to 

A* m = * m A m i m = l,2,...,n, (10) 

where A m are real eigenvalues. 

There is an important difference between the structure of hermitian operators in complex and 
quaternionic quantum mechanics. In complex quantum mechanics we can always trivially relate an 
anti-hermitian operator, A, to an hermitian operator, iJ, by removing a factor i, i.e. A = i H . In gen- 
eral, due to the non-commutative nature of the quaternionic field, this does not apply to quaternionic 
quantum mechanics. 
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Let {'I'm} be a set of normalized eigenvectors of A with complex imaginary eigenvalues {i\ m }. 
The anti-hcrmitian operator A is then represented by 



A = *r Ki , 

r=l 

where ^ + := 'J . It is easy to verify that 

A^ m = ^2 *r A r « * m = ^ * r XJ S rm = * ro A„ 



(11) 



In quaternionic quantum mechanics with quaternionic geometry j36[ , the observable corresponding to 
the anti-hcrmitian operator A is represented by the following hcrmitian quaternionic linear operator 



(12) 



The action of this operators on the eigenvectors $> m gives 

H * ro = * m A m . 

The eigenvalues of the operator H are real and eigenvectors corresponding to different eigenvalues are 
orthogonal. 

How to relate the hermitian operator H to the anti-hermitian operator A? A simple calculation 
shows that the operators LiH and HLi does not satisfy the same eigenvalue equation of A. In fact, 



L, H 



LAJ2^r x r 



and 



*r A r j Lj 



: Y *r A r * r f = t * m A„ 



These problems can be avoided by using the right operator i?; instead of the left operator L;. In fact, 
the operator HR t satisfies the same eigenvalue equation of A, 



HR i q„ 




The eigenvalues of the operator — ARi are real and eigenvectors corresponding to different eigenvalues 
are orthogonal. The right hermiticity of this operator is recovered within a quaternionic formulation 
of quantum mechanics based on complex geometry |m| . 

When the space state is finite-dimensional, it is always possible to form a basis with the eigenvectors 
of the operators H and — ARi. When the space state is infinite-dimensional, this is no longer necessarily 
the case. So, it is useful to introduce a new concept, that of an observable. By definition, the hermitian 
operators H or —ARi are observables if the orthonormal system of vectors forms a basis in the state 
space. 

In quaternionic quantum mechanics with quaternionic geometry, the hermitian operator corre- 
sponding to the anti-hermitian operator A of equation (|ll]) is thus given by the operator H of equa- 
tion (O). By adopting a complex geometry, observables can also be represented by complex linear 
hcrmitian operators obtained by multiplying the corresponding anti-hermitian operator A by — R t . 
We remark that for complex eigenvectors, the operators LiH, HLi, HR t and A reduce to the same 
complex operator 



' H = i X r *J 
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We conclude this section by giving an explicit example of quaternionic hermitian operators in a 
finite two-dimensional space state. Let 

(13) 



3j i _ 

be an anti-hermitian operator. An easy computation shows that the eigenvalues and the eigenvectors 
of this operator are given by 



{2i , 4i} and | ^= 



M JL I k 

j h V2 I 1 



It is immediate to verify that iA and Ai are characterized by complex eigenvalues and so cannot 
represent quaternionic observables. In quaternionic quantum mechanics with quaternionic geometry, 
the quaternionic observable corresponding to the anti-hermitian operator of equation (^) is given by 
the hermitian operator 

F = ^2^ + * 2 4*t= ( 3 fc ^ . (14) 

Within a quaternionic quantum mechanics with complex geometry, a second equivalent definition of 
the quaternionic observable corresponding to the anti-hermitian operator of equation (|l^) is given by 
the complex linear hermitian operator 

-i 3j 



H= ( 3 - ^ ) R, . (15) 

5 THE QUATERNIONIC SCHRODINGER EQUATION 

For simplicity, we shall assume a one-dimensional description. In the standard formulation of quan- 
tum mechanics, the wave function of a particle whose potential energy is V(x,t) must satisfy the 
Schrodinger equation 

ihdMx, t)=H t) = [-£d„ + V(x, t)] t) . (16) 
Let us modify the previous equation by introducing the quaternionic potential 

[V + h-V] (x,t) . 

The i-part of this quaternionic potential violates the norm conservation. In fact, 

The j/fc-part of h ■ V is responsible for T- violation B. To show that, we briefly discuss the time 
reversal invariance in quaternionic quantum mechanics. The quaternionic Schrodinger equation in 
presence of a quaternionic potential which preserves norm conservation, is given by P,p^|,[lq,|36| 

ihd t $(x,t) = [H-jW] $(x,t) , (17) 

where W € C. Evidently, quaternionic conjugation 

-hd t $(x,t)i = H$(x,t) +$(x,t)jW 

does not yield a time-reversed version of the original Schrodinger equation 

-ifrd t <f> T (x,-t) = [H-jW] $ T (x,-t) . (18) 
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To understand why the T-violation is proportional to the j//c-part of the quaternionic potential, let 
us consider a real potential W. Then, the Schrodinger equation has a T-invariance. By multiplying 
the equation ([l7|) by j from the left, we have 

-ihd t j$(x,t) = [H-jW] j$(x,t) , W € M , 

which has the same form of equation (|l8|). Thus, 

$ T (x,-t) =j$(x,t) . 

A similar discussion applies for imaginary complex potential W £ i R. In this case, we find 

$ T (x,-i) = k$(x,t) . 

However, when both V 2 and V 3 are non zero, i.e W £ C, this construction does not work, and the 
quaternionic physics is T-violating. The system of neutral kaons is the natural candidate to study the 
presence of effective quaternionic potentials, V + h ■ V . In studying such a system, we need of V l and 
V 2 , 3 in order to include the decay rates of K s /K^ and CP-violation effects. 



5.1 Quaternionic stationary states 

For stationary states, 

V(x,t)=V(x) and W(x, t) = W{x) , 
we look for solutions of the Schrodinger equation of the form 

= *(*)£(*) • (19) 
Substituting ( |l9| ) in the quaternionic Schrodinger equation, we obtain 

ihV(x) C(f) = [H-j W(x) ] *(x) ((t) . (20) 
Multiplying by — ^ (x) i from the left and by from the right, we find 

ftC(*)?(*)/lC(*)l" = ^) [-iW + fcW(a;)] Ht(x)/\*(x)\* . (21) 

In this equation we have a function of t in the left-hand side and a function of x in the right-hand 
side. The previous equality is only possible if 

ftC(«)C(*)/K(«)| a = *(s) [-*W + fcW(aO] ^/^(x)]^^ (22) 

where g is a quaternionic constant. The energy operator —iTi + k W(x) represents an anti-hermitian 
operator. Consequently, its eigenvalues are purely imaginary quaternions, q = h ■ E. By applying the 
unitary transformation u, 



H E u = -iE , E = ^E\ + El + E. 



3 ' 



equation (22) becomes 

hu((t)((t)u/\((t)\ 2 =u¥(x) [-iH + kW(x)\ ^(x)u/\^{x)\ 2 = -iE . (23) 

The solution $(a;, t) of the Schrodinger equation is not modified by this similarity transformation. In 
fact, 

$(x,t) i&(x)uu((t) = $>{x)((t) . 

By observing that |$(a;,£)| 2 = |^(a;)| 2 |^(t)| 2 , the norm conservation implies |C(0| 2 constant. Without 
loss of generality, we can choose |C(^)| 2 = 1- Consequently, by equating the first and the third term in 
equation (23) and solving the corresponding equation, we find 

C(t) = exp[-iEt/h] C(0) , (24) 

with £(0) unitary quaternion. Note that the position of £(0) in equation ( p4| ) is very important. In 
fact, it can be shown that £(0) exp[— iEt/H] is not solution of equation (^3|). Finally, to complete the 
solution of the quaternionic Schrodinger equation, we must determine ^(x) by solving the following 
second order (right complex linear) differential equation 



2m 



d xx -iV{x) + kW{x) =-^{x)iE . (25) 
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• Real potential 

For W(x) — 0, equation (p5| ) becomes 

£ d„ - V(x) 1 {[9(x)] c - j b'*(z)]c} = i {[*(s)]c - J L?*(*)]c} * E 



Consequently, 



and 



[*( !B )] c = -[*( a; )] c £ 



#d„-V{x)} \j*(x)] c = \j*(x)] c E 



By solving these complex equations, we find 



fc 2 + 



*(x)=exp Jf? {V - E) x fc 1 + exp --y/f? (V-E); 

j { exp [^f? (V+^as] fc 3 + exp [-/§? (y + S ) X J k 4 } 
where n = 1, ...,4, are complex coefficients determined by the initial conditions. 



(26) 



• Free particles 

For free particles, V(x) — W(x) — 0, the previous solution reduces to 

^(x) = exp [i | x] k 1 + exp [—i^x] k 2 + j {exp [fx] k 3 + exp [—fa;] fc 4 } , 



where p = \j2mE. For scattering problems with a wave function incident from the left on quaternionic 
potentials, we have 



^(x) — exp[i j x] + r exp[ 



■jf exp[fx 



(27) 



where \r\ 2 is the standard coefficient of reflection and |f exp[ § a;]| 2 represents an additional evanescent 
probability of reflection. In our study of quaternionic potentials, we shall deal with the rectangular 
potential barrier of width a. In this case, the particle is free for x < 0, where the solution is given by 
(^7j), and x > a, where the solution is 

&(x)=texp[i$x]+jtexp[-%x]. (28) 

Note that, in equations (^7|) and (|2|), we have respectively omitted the complex exponential solution 
exp [—fx] and exp[f x] which are in conflict with the boundary condition that ty(x) remain finite 
as x — > — oo and x — * oo. In equation (|28|), we have also omitted the complex exponential solution 
exp[— i § x] because we are considering a wave incident from the left. 



6 QUATERNIONIC LINEAR DIFFERENTIAL EQUATION 

Consider the second order quaternionic linear differential operator 

T> B = d xx + (a +L-a)d x + b + L-b eA M ®(D. 
We are interested in finding the solution of the quaternionic linear differential equation 

V B <p(x) = . 

In analogy to the complex case, we look for solutions of exponential form 

ip(x) — exp[gx] , 



(29) 



where gel and x € M. To satisfy equation (|29|), the constant q has to be a solution of the quaternionic 
quadratic equation p§-[4l| 



q 2 + (a + h ■ a) q + b a + h ■ b = 



(30) 
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6.1 Quaternionic quadratic equation 

To simplify our discussion, it is convenient to modify equation ( p0| ) by removing the real constant a . 
To do this, we introduce a new quaternionic constant p defined by p — q+ The quadratic equation 
( |30| ) then becomes 

p 2 + ft ■ ap + c + ft ■ c = , (31) 

where c„ = b„ — ^9- and c = b — ^ a. We shall give the solution of equation (|3l]) in terms of real 
constant c„ and of the real vectors a and c. Let us analyze the following cases: 

(i) a x c = 0, 

• a 7^ 0, c 7^ 0: (ii) a • c = 0, 

(iii) axc/O^d'C; 

• a = 0, c ^ 0; 

• a 7^ 0, c = 0; 

• a = c = 0. 

• (i) a x c = 0. In this case a and c are parallel vectors, so equation (|3|) can be easily reduced 
to a complex equation. In fact, by introducing the imaginary unit X = h ■ a/\a\ and observing that 
h ■ c = X a, with a £ K, we find 

p 2 +l\a\p + c + la = , 
whose complex solutions are immediately found. 

• (ii) a ■ c = 0. By observing that a, c and axe are orthogonal vectors, we can rearrange the 
imaginary part of p, h ■ p, in terms of the new basis (a, c, a x c), i.e. 

j) = p G + /i'(ia + i/c+zaxc) . (32) 

Substituting (^) in equation (|3l]) , we obtain the following system of equations for the real variables 
Po, y and z, 

M : pg - (x 2 + x) \a\ 2 - y 2 \c\ 2 - z 2 |a| 2 |c| 2 + c = , 

ft, • a : p Q (1 + 2x) = , 

ft • c : 1 + 2p y — z |a| 2 = , 

ft-axc: y + 2p o 2 = 0. 

The second equation, p (1 + 2x) = 0, implies p = and/or .x = —5. For p n = 0, it can be shown 
that the solution of equation (|3l|), in terms of p , a;, y and z, is given by 



Po = 0, x = -i±VA, y = 0, 2=i, (33) 

a 



where 



For a; = — i, we find 



and 



It is easily verified that 



4 a 2 V u a 



2p„ 



4p 2 + |a| 2 4p 2 



(34) 



v 2 = - 



±2 v /c 2 + |c| 2 -2c -|a| : 



A<0 Vc 2 + |c| 2 - c > i^- 



10 



Stefano De Leo, Gisele Ducati: QUATERNIONIC DIFFERENTIAL OPERATORS 



thus 



a- 



Summarizing, for A ^ 0, we have two quaternionic solutions, p ± =/= p 2 , 



A > 



Po = , 
= -I±VA, 

y = o, 
i . 
W > 



z 



A<0 : p =±i, 2 L^+H 



_ _i 

2 ' 
2po 

1 



4p^ + |a|2 • 

For A = 0, these solutions tend to the same solution p 1 = p 2 given by 



(35) 



(36) 



(37) 



A = 0: p o = 0, x=—i, y = , z = 



(38) 



• (iii) aXc^O/o'C. In discussing this case, we introduce the vector d = c — d a, d = a ■ c/\a\ 2 
and the imaginary part of p in terms of the orthogonal vectors a, d and a x d, 

P = Pa + h-(xa + yd + zaxd) . (39) 

By using this decomposition, from equation ( j3l"l) we obtain the following system of real equations 

K : p 2 - (x 2 + x) |a| 2 - y 2 |d| 2 - z 2 |a| 2 |d| 2 + c = , 

h a : p (l + 2x)+d =0 , 

h<2; 1 + 2p y - z |a| 2 = , 

h ■ a x d : y + 2 p z = . 

The second equation of this system, p (1 + 2x) + d = 0, implies p 7^ since rf ^ 0. Therefore, we 
have 



and 



Po + d 
x = , y 



16w 3 + 8[|a| 2 + 2c ] w 2 +4 



2 Po 



1 



Apl + |a| 2 ' ■ 4p 2 + \a\ 



\a\ 2 (c -dl) + ^-\d\ 2 ,r-d;Aa 1 = 11 



(40) 



(41) 



where w = p 2 . By using the Descartes rule of signs it can be proved that equation ( |4l"| ) has only one 
real positive solution |38j, u> = a 2 , a € R. This implies po = ± a. Thus, we also find two quaternionic 
solutions. 

• a = and c^0. By introducing the imaginary complex unit I = h ■ c/\c\, we can reduce equation 
( |3l| ) to the following complex equation 

p 2 +c +X|c| =0 . 

• a ^ and c = 0. This case is similar to the previous one. We introduce the imaginary complex unit 
X = h ■ a/\a\ and reduce equation (pi]) to the complex equation 



p 2 +I\a\p + c = 



Stefano De Leo, Gisele Ducati: QUATERNIONIC DIFFERENTIAL OPERATORS 



11 



• a = c = 0. Equation (|3l]) becomes 

P 2 + c„ = . 

For c = —a 2 , a e 1, we find two real solutions. For c = a 2 , we obtain an infinite number of 
quaternionic solutions, i.e. p — h ■ p, where \p\ — \a\. 

Let us resume our discussion on quaternionic linear quadratic equation. For a = and/or c = 
and for a x c = we can reduce quaternionic linear quadratic equations to complex equations. For 
non null vectors satisfying a-c = Qoraxc=/=Q=/=a - c, we have effective quaternionic equations. In 
these cases, we always find two quaternionic solutions (|3^), ( p7| ) and (|4C|-^l|). For a ■ c = and A = 0, 
these solutions tend to the same solution ( |38| ) . Finally, the fundamental theorem of algebra is lost for 
a restricted class of quaternionic quadratic linear equations, namely 

q 2 + a 2 = , a € E . 



6.2 Second order quaternionic differential equations with constant coefficients 



Due to the quaternionic linearity from the right of equation (|29|), we look for general solutions which 
are of the form 

tp(x) = <Pi(x) c x + ip 2 (x) c 2 , 

where <Pi{x) and (p 2 (x) represent two linear independent solutions of equation ( |29| ) and c x and c 2 are 
quaternionic constants fixed by the initial conditions. In analogy to the complex case, we can distin- 
guish between quaternionic linear dependent and independent solutions by constructing a Wronskian 
functional. To do this, we need to define a quaternionic determinant. Due to the non-commutative na- 
ture of quaternions, the standard definition of determinant must be revised. The study of quaternionic, 
complex and real functionals, extending the complex determinant to quaternionic matrices, has been 
extensively developed in quaternionic linear algebra |^6|-^1 . In a recent paper |5Q] , we find an interest- 
ing discussion on the impossibility to obtain a quaternionic functional with the main properties of the 
complex determinant. For quaternionic matrices, M, a real positive functional, |detM| = y^det [MM^], 
which reduces to the absolute value of the standard determinant for complex matrices, was introduced 
by Study [5l| and its properties axiomatized by Dieudonne |p2[ . The details can be found in the excel- 
lent survey paper of Aslaksen H] . This functional allows to construct a real positive Wronskian [fjll 



dot 



Solutions of equation 



tpi(x) <p 2 {x) 
<Pi(x) <p 2 (x) ^ 

<Pi{x)\ \<£ 2 {x) - ip^x) <p~ x {x) p 2 (x)\ 
ip 2 (x)\ \ipiix) - <p 3 (x) tp^ix) tp x {x)\ 
<P-l{x)\ \<p 3 (x) - (p^x) tp~ x {x) <p 2 (x)\ 

(f 2 (x)\ \<P!(X) - Lp 2 (x)lp 2 1 (x)Lp 1 (x)\ . 



ip ll2 (x) = exp[q lt2 x] = exp[(p 1>2 ~ 9 f)x] 



are given in terms of the solutions of the quadratic equation (|3l|), p 1)2 , and of the real variable x. In 
this case, the Wronskian becomes 

W(x) = \Pi - Pal lexp^x]! |exp[g 2 x]| . 

This functional allows to distinguish between quaternionic linear dependent (W = 0) and independent 
(W 7^ 0) solutions. A generalization for quaternionic second order differential equations with non 
constant coefficients should be investigated. 

For p t ^ p 2 , the solution of equation ( p9|) is then given by 



ip(x) = exp[— ^ x] {exp[p 1 x] c t + exp[p 2 x] c 2 } 



(42) 



12 



Stefano De Leo, Gisele Ducati: QUATERNIONIC DIFFERENTIAL OPERATORS 



As observed at the end of the previous subsection, the fundamental theorem of algebra is lost for a 
restricted class of quaternionic quadratic equation, i.e. p 2 + a 2 — where aei For these equations 
we find an infinite number of solutions, p = h ■ a with \a\ 2 — a 2 . Nevertheless, the general solution 
of the second order differential equation 



<p(x) + a 2 <p(x) = , 
is also expressed in terms of two linearly independent exponential solutions 

tp(x) = exp[z ax]c ± + exp[— i ax]c 2 ■ 



(43) 
(44) 



Note that any other exponential solution, exp[h-a x], can be written as linear combination of exp[« ax] 
and exp[— iax], 

exp[fa • a x] = y~ {exp[i ax] (a ~ i h ■ a) + exp[— i ax] (a + i h ■ a)} . 

As consequence, the loss of the fundamental theorem of algebra for quaternions does not represent an 
obstacle in solving second order quaternionic linear differential equations with constant coefficients. 
To complete our discussion, we have to examine the case p x — p 2 ■ From equation (|38|) we find 



Pi = Vi 



-tKk h ■ a x (b 

\a\ 2 \ 



Thus, a first solution of the differential equation (29) is 



^)=cxp{[fc.(^-f)-f] x} 



For a x b = 0, we can immediately obtain a second linearly independent solution by multiplying 
exp[— f x] by x, r/(x) = x£(x). For a x b ^ 0, the second linearly independent solution takes a more 
complicated form, i.e. 



r)(x) 



ha 

W 



(45) 



It can easily be shown that rj(x) is solution of the differential equation (p9[), 

rj(x) + a r)(x) + b T)(x) = x (q 2 + aq + b) + 2 q + a + ^ {q 2 + a q) + 



bM 
m 



2q + a- 
(2h- 2$. 



b, 



ha 

W 



h b, 

' \a\ 2 



= . 



Thus, for p ± = p 2 = p = h 

by 



— j ) , the general solution of the differential equation (|29|) is given 



ip(x) = exp[— t& x] |cxp[pa;] c x 



h-a 

w 



exp 



[px] c 2 | 



(46) 



6.3 Diagonalization and Jordan form 



To find the general solution of linear differential equations, we can also use quaternionic formulations 
of eigenvalue equations, matrix diagonalization and Jordan form. The quaternionic linear differential 
equation (E9h can be written in matrix form as follows 



(47) 



where 



M = 



1 



and $(x) 



tp(x) 
(p(x) 
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By observing that x is real, the formal solution of the matrix equation (^) is given by 

$(x) = exp[Af x] $(0) , 



(48) 



where $(0) represents a constant quaternionic column vector determined by the initial conditions 

" , (Mi) 



ip(0), 0(0) and exp[M x] = Y^ - 



In the sequel, we shall use right eigenvalue equations for 
quaternionic linear matrix operators equations 



(49) 



Without loss of generality, we can work with complex eigenvalue equations. By setting 'J = $m, from 
the previous equation, we have 



Mt = M$i/ = $gu = §u uqu = * z , 



(50) 



where z € C and u is a unitary quaternion. In a recent paper ]3l| ] , we find a complete discussion of the 
eigenvalue equation for quaternionic matrix operators. In such a paper was shown that the complex 
counterpart of the matrix M has an eigenvalue spectrum characterized by eigenvalues which appear 
in conjugate pairs {z 1 ,z x ,z 2 ,Z 2 }. Let ^ and ^ 2 be the quaternionic eigenvectors corresponding to 
the complex eigenvalues z t and z 2 

M *! = Zi and M * 2 = * 2 z 2 . 

It can be shown that for \z t \ ^ \z 2 \, the eigenvectors and \& 2 are linearly independent on H and 
consequently there exists a 2 x 2 quaternionic matrix S = [ i Sf 1 ^ 2 ] which diagonalizes M, 



exp [Ma;] = S exp 



z t 
z 2 



s- 1 = s 



( exp^x] 
I exp[2 2 x] 



s- 1 



In this case, the general solution of the quaternionic differential equation can be written in terms of 
the elements of the matrices S and S~ x and of the complex eigenvalues z x and z 2 , 



<p(x) 



S 1± exp^cc] S 12 exp[z 2 x] 
S 2 i exp^x] exp[z 2 a;] 



S-V(0) + 5-^(0) 



Hence, 



(p(x) 



S u expfc x] [S- 1 p(0) + S- 1 0(0)] + 
S 12 exp[z 2 x] [S-V(0) + -5-^(0)] 
expfS^ (Su)" 1 1] S„ [5 n >(0) + 5" 1 0(0)] 
exp [S„ z 2 (S ia ) -1 S M p(0) + S- 1 0(0)] 
exp [S„ (S^)- 1 x] Sn [S- 1 <p(0) + S- 1 0(0)] + 
exp [S 22 (S^y 1 x] S 12 [S- 1 <p(0) + S- 1 0(0)] . 



(51) 



We remark that a different choice of the eigenvalue spectrum does not modify the solution (pl|) . In 
fact, by taking the following quaternionic eigenvalue spectrum 



{qi , q 2 } = {u 1 z 1 u 1 , u 2 z 2 u 2 } , \q x \ ^ \q 2 \ , 
and observing that the corresponding linearly independent eigenvectors are given by 

{$ 1 = $ 2 = * 2 w 2 } , 

we obtain 

A/ = [$ t $ 2 ] diag {?!,&} [*x $ a ] -1 

= $ 2 « 2 ] diag{uiZiMi,S 2 ^ a iij,} [vfiUi * 2 it 2 ] _I 

= diag{z 1 ,z 2 } 4-,]- 1 . 



(52) 
(53) 
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Let us now discuss the case \z t \ = \z 2 \. If the eigenvectors {$„, ^ b }, corresponding to the eigenvalue 
spectrum { z , z }, are linearly independent on H, we can obviously repeat the previous discussion and 
diagonalize the matrix operator M by the 2x2 quaternionic matrix U — [*i> 1 ^ 2 ]■ Then, we find the 

<p(x) = exp [U u z (U^y 1 x] U n [U- 1 <p(Q) + U£ 0(0)] + 
exp [U 12 z {U^Y 1 x] U 12 [U-, 1 tp(0) + U£ 0(0)] 
= exp [U 21 (U^y 1 x] U u [U- 1 <p(0) + U- 1 0(0)] + 

exp [U 22 (U^y 1 x] U» [U 2 y <p(0) + U- 1 0(0)] . (54) 

For linearly dependent eigenvectors, we cannot construct a matrix which diagonalizes the matrix 
operator M. Nevertheless, we can transform the matrix operator M to Jordan form 



m=j(i 

It follows that the solution of our quaternionic differential equation can be written as 



(55) 



<5>(x) — J exp 



z 1 
z 



J 11 X J±i ~\~ J\2 



J- x $(0) 
exp[za;] 



j- 1 ?(o) + Jr, 1 0(0) 

J- 1 <p(0) + J- 1 0(0) 



Thus, 



<p(x) = J„ exp[zx] [J- 1 <p(p) + J" 1 0(0)] + 

(a:J u + J 12 ) expire] [J^ 1 <p(0) + J^ 1 0(0)] 
= exp [J u z {J,,)- 1 x] J u [J- 1 ip(0) + J" 1 0(0)] + 

[a; + J 12 (Jii) -1 ] exp[JnZ (Jn) _1 a;] x 
J„ [J,? V (0) + J,- 1 0(0)] 
= exp [J 21 {J^y 1 x] J u [J" 1 p(Q) + J" 1 0(0)] + 
[x + J 12 (Jn) -1 ] exp [ J 21 (Ju) -1 x] x 
J„ [J,- 1 ^(0) + J,- 1 0(0)] . 



(56) 



The last equality in the previous equation follows from the use of equation (^5|) and the definition of 
M. 

Finally, the general solution of the quaternionic differential equation (|2^ ) can be given by solv- 
ing the corresponding eigenvalue problem. We conclude this section, by observing that the quater- 
nionic exponential solution, exp[ga;], can also be written in terms of complex exponential solutions, 
u exp[zx\ u~ x , where q = uzu" 1 . The elements of the similarity transformations S, U or J and the 
complex eigenvalue spectrum of M determine the quaternion u and the complex number z. This 
form for exponential solutions will be very useful in solving complex linear differential equations with 
constant coefficients. In fact, due to the presence of the right acting operator we cannot use 
quaternionic exponential solutions for complex linear differential equations. 



7 COMPLEX LINEAR QUATERNIONIC DIFFERENTIAL EQUATIONS 

Consider now the second order complex linear quaternionic differential operator 

T> c = [a 02 + L-a 2 + (b 02 + L ■ b 2 ) R,] d m + 
[a 01 +L-a, + (b 01 +L-b 1 ) R t ] <9 X + 
a 00 + L ■ a + (6 00 + L ■ b ) R % 
G Ac ® O 
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and look for solutions of the complex linear quaternionic differential equation 

V c <p(x) = . 



(57) 



Due to the presence of Ri in ([57]), the general solution of the complex linear quaternionic differential 
equation cannot be given in terms of quaternionic exponentials. In matrix form, equation (p7p reads 



$(a>) = M c <S>(x) 



(58) 



where 



1 

■b c -a r 



and 



<p(x) 



The complex counterpart of complex linear quaternionic matrix operator M c has an eigenvalue spec- 
trum characterized by four complex eigenvalues {z u z 2 , z 3 , z 4 }. It can be shown that M c is diagonal- 
izable if and only if its complex counterpart is diagonalizable. For diagonalizable matrix operator 
M c , we can find a complex linear quaternionic linear similarity transformation S c which reduces the 
matrix operator M c to diagonal form plj 



M c = S c 
It is immediate to verify that 



2 







Z3+Z4 

2 



are eigenvectors of the diagonal matrix operator 



2 



^^R % 







2 1 



Z3+Z4. 
2 



Z3 — Z4 D 

2% Li 



Z3 — Z4. JD 

2i 1Li 



with right complex eigenvalues z ± , z 2 , z 3 and z 4 . The general solution of the differential equation (57) 
can be given in terms of these complex eigenvalues, 



<p(x)=S cl t exp[(^ + 



^Ri) x\ [5- 1 1 1 ^(0) + ^: 1 1 2 ^(0)] 

2i 



S C12 exp [(^ + R<) x] ^(0) + S-> 2 #0)] 
«i exp[zi x] kx + u 2 exp[z 2 a;] fc 2 + 
u 3 exp[z 3 x] k 3 + u 4 exp[z 4 x] fc 4 , 



(59) 



where k n are complex coefficients determined by the initial conditions. This solution holds for diago- 
nalizable matrix operator M c . For non diagonalizable matrix operators we need to find the similarity 
transformation J c which reduces M c to the Jordan form. For instance.it can be shown that for equal 
eigenvalues, z 1 = z 2 , the general solution of the differential equation ( |57|) is 



<p(x) = u exp[z x] k t + (u x + u) exp[z x] k 2 + u 3 exp[z 3 x] k 3 + u 4 exp[z 4 x] k 4 



(60) 



7.1 Schrodinger equation 

Let us now examine the complex linear Schrodinger equation in presence of a constant quaternionic 
potential, 



£d„-V + jW 9(x) = iV(x)iE 



(61) 



In this case, the complex linear matrix operator 

1 



M c = 



■be 



, b c = V — jW + iERi , 
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represents a diagonalizable operator. Consequently, the general solution of the Schrodinger equation 
is given by 



tp(x) = u 1 exp[z! x] k 1 + u 2 exp[z 2 x] k 2 + u 3 exp[z 3 x] k 3 + it 4 exp[z 4 x) fc 4 . 



(62) 



The quaternions u n and the complex eigenvalues z n are obtained by solving the eigenvalue equation 
for the complex linear operator M c . We can also obtain the general solution of equation d6~T|) by 



substituting u exp^ ^ zx] in the Schrodinger equation. We find the following quaternionic equation 

uz 2 ~(V-jW)u-iEui = , 
where u = z u + j z u . This equation can be written as two complex equations 

[z 2 -(V- E)] z u - Wz u = [z 2 -(V + E)] z u + W z u = Q . 
An easy calculation shows that z satisfies the complex equation 

z i -2Vz 2 + V 2 + \W\ 2 -E 2 = 0, (63) 

whose roots are 



z M =±y/V- V E2 - \W\ 2 =±z_ and z 3 , 4 = ± yJV + V ' E 2 - \W\ 2 = ± z. 
By setting (ti 1;2 ) c = (-ju St4 ) c = 1, we find 

and u + 



1+3 ■ 



E+^E 2 -\W\ 2 



+ 3 



-B+V 'E 2 -\W\' 

The solution of the complex linear quaternionic Schrodinger equation is then given by 

\&(x) = u_ |exp[^/^z_ x] ki + exp[- ^f? 2- x] k 2 \ + 
u + |exp[^/^z + x] k 3 + exp[- \J~^ z+ x]kA . 



(64) 



(65) 



(66) 



Equation (63) can also be obtained by multiplying the complex linear Schrodinger equation (|6lj) from 
the left by the operator 



$- d mx - V - j W 

2m xx J 



This gives 

V 2m J xxxx 



-2 h 



2 m 



Vd^+V 2 + \W\ 2 *(x) 



2m x ' 

E 2 ^(x) 



V + jW 



^{x)iE 



By substituting the exponential solution u exp[y'^za;] in the previous equation, we immediately 
re-obtain equation (|63]). 



8 QUATERNIONIC CONSTANT POTENTIALS 

Of all Schrodinger equations the one for a constant potential is mathematically the simplest. The 
reason for resuming the study of the Schrodinger equation with such a potential is that the qualitative 
features of a physical potential can often be approximated reasonably well by a potential which is 
pieced together from a number of constant portions. 



Stefano De Leo, Gisele Ducati: QUATERNIONIC DIFFERENTIAL OPERATORS 



17 



8.1 The potential step 

Let us consider the quaternionic potential step, 

V(x) - j W(x) = 



x < 
V - j W x>Q 



where V and W represent constant potentials. For scattering problems with a wave function incident 
from the left on the quaternionic potential step, the complex linear quaternionic Schrodinger equation 
has solution 



( x < : 

exp[?|x] 



r exp[— i ? x] 



■jf exp[|x] 



x > : 

u_ t expL/ z_ x] + u + t exp[- 



(67) 



u_ t exp[- 



z_ + u + t exp[- 



/ 2 m 



-2+ 



[E < ^JV 2 + \W\ 2 } 



where r, f, t and t are complex coefficients to be determined by matching the wave function and 
its slope at the discontinuity of the potential x — 0. 

For E > \/V 2 + \W\ 2 , the complex exponential solutions of the quaternionic Schrodinger equation 
are characterized by 



and 



+ = ^E 2 -\W\ 2 + V 



W z + X ] 



z-=i\j sjE 2 - \W\ 2 -V e i I 

The complex linearly independent solutions 

u_ exp[— j£ z_ x] and u + exp[ 

have been omitted, k 2 = k 3 = in (66), because we are considering a wave incident from the left 

and because the second complex exponential solution, exp[y / '~^J z + a;], is in conflict with the boundary 

condition that \l/(x) remain finite as x — > oo. The standard result of complex quantum mechanics 
are immediately recovered by considering W — and taking the complex part of the quaternionic 

solution. 

For E < \/V 2 + \W\ 2 , the complex exponential solutions of the quaternionic Schrodinger equation 
are characterized by 



z_ = ^ V -^E 2 -\W\ 2 , 
z ± = (V 2 + \W\ 2 - E 2 )i exp[±if] , tan 6 

The complex linearly independent solutions 



tJ\w\ 2 -e 2 
v 



6 



[E > \W\] 
[E<\W\] 



u_ exp[y' z_ x] and u + exp[y z + x] 

have been omitted, k t — k 3 — in (p6|), because they are in conflict with the boundary condition that 
^(x) remain finite as x — > oo. 

A relation between the complex coefficients of reflection and transmission can immediately be 
obtained by the continuity equation 



where 



and 



d t p(x,t) + d x J(x,t) = , 
p(x,t) =$(x,t)$(x,t) , 
t) = -t{ [^M^ t)] i *(x, t) - $(x, t) i d.Q(x, t) } 
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Note that, due to the non commutative nature of the quaternionic wave functions, the position of 
the imaginary unit i in the probability current density J(x,t) is important to recover a continuity 
equation in quaternionic quantum mechanics. For stationary states, = \I/(x)exp[ — i 7^ i]C(0), 

it can easily be shown that the probability current density 



J(x,t) = ±C(0) exp[if t}{ [0.*(a;)] i*(x)-*(x)id x y(x)} exp[-if t]C(0) . 
must be independent of x, J(x,t) — f(t). Hence, 

£ { [&*(*)] i*(x)-*(x)id.*(x)} =exp[-if *]C(0)/(i)C(0) exp[if t] =a. 
where a is a real constant. This implies that the quantity 

J = & { [dMx)] * - *(*) i } 
has the same value at all points x. In the free potential region, x < 0, we find 

^- = £ (i-H 2 ) • 

In the potential region, x > 0, we obtain 



J + 



>± ( j e * - \w\ 2 - v) \i-( — V 



\t\ 2 [E> y/V* + \W\>] 
[E< ^/V 2 + \W\ 2 ] 



Finally, for stationary states, the continuity equation leads to 

1 



y/E 2 -\W\ 2 -V 

E 



\W\ 



E+y/E 2 -\W\' 



\t\ 2 = l [E> y/V* + \W\ 3 ] 



(69) 



1 



[E<y/V' + \W\'\ 



Thus, by using the concept of a probability current, we can define the following coefficients of trans- 
mission and reflection 



R= \r\ 2 , T 
R= \r\ 2 , T = 



y/E 2 -\W\ 2 -V 

E 



1 



\W\ 



E+,/E 2 -\W\- 



\t\ 2 [e> jw+WY 2 ] , 



[E< y/V' + \W\*] . 



These coefficients give the probability for the particle, arriving from x — — oo, to pass the potential 
step at x = or to turn back. The coefficients R and T depend only on the ratios E/V and |W|/V". 
The predictions of complex quantum mechanics are recovered by setting W = 0. 



8.2 The rectangular potential barrier 

In our study of quaternionic potentials, we now reach the rectangular potential barrier, 



V(x) - j W(x) 



x < 
V -jW <x<a 
x > a 



For scattering problems with a wave function incident from the left on the quaternionic potential 
barrier, the complex linear quaternionic Schrodinger equation has solution 

r x < : 

exp[i | x] + r exp[— i \ x\ + jf exp[ fx]; 



< x < a 



*(x) = < 



Jcxpiyi? z_x] k.+expl- tJ%£z_x] k 2 \ + 
expL/fj z + x] k 3 + exp[- 



(70) 



T5T Z + X] k 4 



x > a : 

k t exp[i I x] + jt exp[- § x] . 
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The complex coefficients r, f, t and t are determined by matching the wave function ^>(x) and its 
slope at the discontinuity of the potential x — and will depend on \ W\. 

By using the continuity equation, we immediately find the following relation between the trans- 
mission, T — \t\ 2 , and reflection, R = \r\ 2 , coefficients 

R + T = l. (71) 



8.3 The rectangular potential well 

Finally, we briefly discuss the quaternionic rectangular potential well, 



V(x) - j W{x) 



x < 

-V + jW 0<x<a 
x > a 



In the potential region, the solution of the complex linear quaternionic Schrodinger equation is then 
given by 



ty(x) — u_ |exp[y /r 2jj z„ x] ki + exp[- 
u + |exp[^/ z + x] k 3 + exp[- 



x] fc 2 } + 

■}■ 



(72) 



where 



and 



U- = 1- j 



E + y/E 2 -\W\ 2 



E + y/E 2 -\W\ 2 



\W\ 2 + V 



z + = ^E 2 -\W\ 2 ~ . 



Depending on whether the energy is positive or negative, we distinguish two separate cases. If E > 0, 
the particle is unconfined and is scattered by the potential; if E < 0, it is confined and in a bound 
state. We limit ourselves to discussing the case E < 0. For \W\ < \E\ < \JV 2 + \W\ 2 , solution (72) 
becomes 



|exp[i y/fj yJ^/EP - \Wf + Vx] k ± + exp[- i yf^/% 2 - \w\ 2 + V x] fc 2 | + 
|exp[i \p?$-\Jv - sfW- \W\ 2 x\ fc 3 +exp[-i y^V- ^E 2 - \W\ 2 x] kA . 



(73) 



For \E\ < \W\, the solution is given by 



u. 



Jexpl^^ exp[z ttE] x ] kl + ex p[- p exp[z ^] x] fc 2 } + 
eM^~P ex P[* x \ fc 3 + exp[- exp[~ i 2±£] x] fc 4 } , 



(74) 



a/| W\ 2 — 

where p = ^V 2 + \W\ 2 ~ E 2 and tan0 = v ' [ v . In the region of zero potential, by using the 
boundary conditions at large distances, we find 

( x < : 



exp[ 



x > a 



2m 
li 2 



E\x]c 1 + j exp[ — i 



f£\E\x\c 4 ; 



(75) 



exp[-yJ^\E\x]d 2 +j exp[i^f? |£|x]d 3 . 
The matching conditions at the discontinuities of the potential yield the energy eigenvalues. 
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9 CONCLUSIONS 

In this paper, we have discussed the resolution of quaternionic, T> m ip(x) = 0, and complex, T> c (p(x) = 0, 
linear differential equations with constant coefficients within a quaternionic formulation of quantum 
mechanics. We emphasize that the only quaternionic quadratic equation involved in the study of 
second order linear differential equations with constant coefficients is given by equation ( |30[ ) following 
from T> B (p{x) = 0. Due to the right action of the factor i in complex linear differential equations, we 
cannot factorize a quaternionic exponential and consequently we are not able to obtain a quaternionic 
quadratic equation from T> c tp(x) = 0. Complex linear differential equations can be solved by searching 
for quaternionic solutions of the form q exp[zx], where q S H and z £ C. The complex exponential 
factorization gives a complex quartic equation. A similar discussion can be extended to real linear 
differential equations, T> x cp(x) = 0. In this case, the presence of left/right operators L and R in T> F 
requires quaternionic solutions of the form q exp[Ax], where q S H and A g R. A detailed discussion 
of real linear differential equations deserves a further investigation. 

The use of quaternionic mathematical structures in solving the complex linear Schrodinger equa- 
tion could represent an important direction for the search of new physics. The open question whether 
quaternions could play a significant role in quantum mechanics is strictly related to the whole un- 
derstanding of resolutions of quaternionic differential equations and eigenvalue problems. The inves- 
tigation presented in this work is only a first step towards a whole theory of quaternionic differential, 
integral and functional equations. Obviously, due to the great variety of problems in using a non- 
commutative field, it is very difficult to define the precise limit of the subject. 
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APPENDIX A: Quaternionic linear quadratic equations 

In this appendix, we give some examples of quaternionic linear quadratic equations, see cases (i)-(iii) 
of subsection fO, and find their solutions. 

• (i): p 2 + V2{i+j)p-l-2V2{i+j) =0. 

In solving such an equation we observe that a = (\/2, \/2, 0) and c = — (2^,2^/2,0) are parallel 
vectors, c = —2 a. Consequently, by introducing the complex imaginary unit X = (i + j)/y2, we can 
reduce the quadratic quaternionic equation to the following complex equation, 

p 2 + 21p- 1-42 = , 

whose solutions are p 12 = —I ± 2 \/X. It follows that the quaternionic solutions are 



Pi 



±\/2- (l T V2 



i + j 
V2 



• (ii): p 2 +ip + ±k = 0, A = 0. 

We note that a = (1,0,0) and c = (0,0,^) are orthogonal vectors and A = 0. So, we find two 
coincident quaternionic solutions given by 



i h ■ a + h ■ a 



i+j 



• (ii): p 2 + jp + 1 - k = 0, A > 0. 
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In this case, a = (0, 1, 0) and c = (0, 0, —1) are orthogonal vectors, c„ = 1 and A = 1/4. So, 

p o =0, x = -\±\, y = 0, z = l. 

By observing that 

h ■ a = j , hc~ — fc , h ■ a x c = —i , 
we find the following quaternionic solutions 

Pl = -i and p 2 = + ■ 

• (ii): p 2 + kp + j = 0, A < 0. 

We have a = (0, 0, 1), c = (0, 1, 0) and c = 0. Then a ■ c = and A = -3/4. So, 
In this case, 

h ■ a = k , h ■ c = j , h ■ a x c = —i , 
thus, the solutions are given by 

p 1)2 = § (±l-i=Fj-A;) . 

• (Hi): p 2 + ip+l + i + k = 0. 

We have a = (1, 0, 0), c = (1, 0, 1) and c = 1. In this case a • c / 0, so we introduce the quaternion 
d„ + h ■ d = 1 + fc, whose vectorial part d — c d a a = (0, 0, 1) is orthogonal to a. The imaginary part 
of our solution will be given in terms of the imaginary quaternions 

h ■ a — i , h- d — k , h ■ a x d = —j . 

The real part of p is determined by solving the equation 

16^ + 24^-3^-1 = 0. 

The real positive solution is given by pl = j. Consequently, 

Pa = ' X = — 2 T 1 ? J/~"F2' Z = 2 ' 

The quaternionic solutions are 

Pi = 2 (1 — 3i — j — A;) and p 2 = -\ (1 - i + j - k) . 

APPENDIX B: Quaternionic linear differential equations 

We solve quaternionic linear differential equations whose characteristic equations are given by the 
examples (i)-(iii) in the previous appendix. 

• (i): l p(x) + V2(i+j)ip(x)- [l + 2v^(*+j)] <p(x)=0, <p(0) = i, 0(0) = ^. 

The exponential exp[px] is solution of the previous differential equation if and only if the quaternion 
p satisfies the following quadratic equation 

p 2 + V2(i+j)p-l-2V2(i+j) = , 
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whose solutions are given by 



Pl . 2 = ±V2- (lT V2) 



i + j 
V2 



Consequently, 

ip(x) = cxp • 



l y/2 - (l - \/2) x| Cl +exp| - y/2 - (l + \/2~) x J c. 



By using the initial conditions, we find 



<p(ar) = cxp 



i+j 
V2 



cosh 



V2 + 



i+j 
V2 



• (ii): <p(x) + (l + i)<p(x) + ?±^<p(x)=0 , <p(0)=0, 0(0)=- 

We look for exponential solutions of the form <p(x) = cxp[qx] = exp[ (p — ^)x]. The quaternion p 
must satisfy the quadratic equation 



p 2 +ip+ \ k = 



This equation implies 



Thus, 



Pi = I>2 



(fii(x) = cxp 



i+j 



l + i+ j 



The second linearly independent solution is given by 



V2 (x) = ( x + i ) exp 



l + i + j 



By using the initial conditions, we find 

<p(x) = {cxp[qx} + ( x + i ) exp[qx] i } [ 1 + q^ 1 ( 1 + i q ) i ] _1 , 
where g = — ( 1 + i + j ) / 2. 

• (ii): 0(ar) + (2 + i)0( a; ) + (2 + i-fc) ^(z) = , y>(0) = ±=i , 0(0)=.?. 
The exponential solution (^(a;) = exp[gx] = exp[ (p — l)x] leads to 

p 2 + jp+ 1 - k = , 

whose solutions are 

Pi = -i and p 2 = -(i+j) ■ 

Consequently, 

ip(x) — cxp[ — x ] { cxp[ — i x } c t + exp[ — ( i + j ) x ] c 2 } . 
The initial conditions yield 

<p(ar) = exp[- x] { exp[- ix] 3 ~' 2 ~ 2j + exp[- (i + j ) x] ( j - 1 ) } . 
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• (ii): (p{x) + k ip(x) + j ip(x) = , ip(0) = i + k , 0(0) = 1 . 
The characteristic equation is 

p 2 + kp + j = , 

whose solutions are 

p 1<2 = \ {±l-iTj-k) . 
Thus, the general solution of our differential equation reads 

ip(x) = cxp[ 1 ~V'~ fc x \ c i +cxp[- x ]c 2 . 

By using the initial conditions, we obtain 

<p(x) = { exp[ i-'-J-fc x ] + exp[- l±^±t± x] j i±fc . 

• (iii): ^5(a;) + (i-2) y>(a;) + (2 + fc) ¥>(a;)=0, p(0) = , y>(0)=j. 

By substituting <p{x) = cx^[qx] = cxp[ (p + 1 )x] in the previous differential equation, we find 

p 2 + ip+l + i + k = 0. 
The solutions of this quadratic quaternionic equation are 

Pi = | (1 - 3i - j - fc) and p 2 = -\ (1 - i + j - k) . 
So, the general solution of the differential equation is 

ip(x) = exp[ 1 ~ 3i -i- k x ] Cl +cxp[- 1 - i +J- k x]c 2 . 
By using the initial conditions, we obtain 

<p(x) = { cxp[ l^t± x] _ cxp [_ k±LL± X ] } t^H . 

APPENDIX C: Diagonalization and Jordan form 

In this appendix, we find the solution of quaternionic and complex linear differential equations by 
using diagonalization and Jordan form. 

Quaternionic linear differential equation 

By using the discussion about quaternionic quadratic equation, it can immediately be shown that the 
solution of the following second order equation 

(f(x) + (k-i) if(x) - j ip(x) = , 

with initial conditions 

<^(0) = |, 0(O) = l + §, 

is given by 

ip(x) = [x+ | ) cxp[i a:] . 
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Let us solve this differential equation by using its matrix form ( p7| ) , with 

M ~ ( 3 i-k 

This quaternionic matrix can be reduced to its Jordan form 



U = ./( j }]./ 



by the matrix transformation 



J= U i + lJ ' J ~ 1 = l-4r ¥ 

The solution of the quaternionic linear quaternionic differential equation is then given by 

tp(x) = J u exp [iz] [Jf/ (f(0) + J^ 1 0(0)] + 

( x J u + J 12 ) exp [ i x } [J-, 1 (p(0) + J" 1 0(0)] 
— (xJ ±1 + Ji 2 ) exp [ix] 
= (x+ | ) exp[jcc] . 

Complex linear differential equations 

Let us now consider the complex linear quaternionic differential equation 

tp(x) - j fix) i = , 

with initial conditions 

tp(0) = j , 0(0) = k . 
To find particular solutions, we set ip(x) — q exp[zx]. Consequently, 

qz 2 - j qi = . 

The solution of the complex linear second order differential equation is 

tp(x) = 5 [(i + j) exp[— ix] + (j — i) coshx + (k — 1) sinhx] . 
This solution can also be obtained by using the matrix 



1 

jR, 



and its diagonal form 
where 

S c = 

and 



1/ = S ( °. ) S 



c > 



2 T 


2 


1+t+J— * 


l+i-j+fc 


2 


2 J 




1+i—j — k 




2 


1 — i+j — k 1 


l-i-j' + fc 



4 I 1 — f + J — fc |_ 1- 1 j + k rj l + i+j + k 1 + i — j — k r> 

\ 2 ' 2 ' 2 2 * 

The solution of the complex linear quaternionic differential equation is then given by 

<p(x) =S cll exp[-iR iX ] [S-l, <p(0) + S-l 2 0(0)] + 
S C12 exp[ix] [S- 2 \ ^(0)+S , - 2 1 2 0(0)] 
= i {(! - i + 3 ~ k ) exp[-x] - (1 + i - j - k) exp[x]} + 
^jr exp[— ix] . 
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